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Abstract 

We discuss a geometrical interpretation of the Z-invariant Ising model in terms of isoradial 
embeddings of planar lattices. The Z-invariant Ising model can be defined on an arbitrary planar 
lattice if and only if certain paths on the lattice edges do not intersect each other more than once 
or self- intersect. This topological constraint is equivalent to the existence of isoradial embeddings 
of the lattice. Such embeddings are characterized by angles which can be related to the model 
coupling constants in the spirit of Baxter's geometrical solution. The Ising model on isoradial 
embeddings studied recently by several authors in the context of discrete holomorphy corresponds 
to the critical point of this particular Z-invariant Ising model. 
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I. INTRODUCTION 



The most general setting in which the two dimensional Ising model is exactly solvable is 
obtained when the model is invariant under star-triangle deformations of the lattice struc- 
ture. This model was introduced by Baxter 1 ' 2 and is known in the literature as the Z- 
invariant Ising model. 

Recently the Ising model on a certain class of planar lattices, those with isoradial embed- 
dings as defined bellow, has attracted attention in connection with the concept of discrete 
holomorphy 3 ' 4 ' 5 . Mercat 3 conjectured a criticality condition and Kenyon 6 showed that the 
restricted model is solvable. In this paper we will prove the criticality conjecture and bring 
these results into perspective by showing that the Ising model considered by these authors is 
the restriction to the critical point of a Z-invariant Ising model. We will show that isoradial 
embeddings provide a geometrical interpretation of the Z-invariant Ising model which gen- 
eralizes the model's geometrical solution 2 in terms of angles characterizing straight rapidity 
lines. 

In the seminal work of Baxter 2 the Z-invariant Ising model is defined on planar graphs or 
'lattices' G, which are defined using an associated rapidity lattice Q s of planar intersecting 
straight lines such that no three lines intersect at a point. These lattices Q s can be colored 
with two colors in such a way that faces with different colors do not share edges. The Z- 
invariant Ising model is defined on the lattice G with vertices on one of the colored subset 
of faces and edges connecting neighboring same color faces which share a vertex, for an 
example see of Fig. 1. The model coupling constants are determined by parameters, called 
the rapidities, associated with the straight lines in Q s . 

In this paper we will address the converse problem. Namely, given an arbitrary planar 
lattice G, defined without reference to a rapidity lattice, is it possible to assign coupling 
constants to its edges in such a away that the Ising model on G is Z-invariant? This 
problem has been considered previously in the context of quasi-crystals 7,8 , by using either 
the de Bruijn grid method 9 or problem specific rapidity line choices 10 . For general planar 
lattices we will show in Section III that a Z-invariant Ising model can be defined if and only 
if the lattice satisfies a certain topological property: that zig-zag paths on the lattice edges, 
this is paths which alternately turn maximally left and right (see Fig. 2), do not intersect 
each other more than once or self- intersect. If this condition is satisfied a Z-invariant Ising 
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model can be denned using rapidity lines in one to one correspondence with the zig-zag 
paths. 

Zig-zag paths are closely related to de Bruijn's 9 concept of skeleton of a parallelogram 
tilling and have recently been used by Kenyon and Schlenker 11 to discuss the existence of 
isoradial embeddings of planar lattices. A strictly convex isoradial embedding of a planar 
lattice is a particularly symmetric drawing of the lattice such that every vertex is at unit 
distance from the center of the faces to which it belongs. For instance Fig. 3 would be an 
isoradial embedding of the lattice shown if drawn in such a way that all edges in dotted line 
have equal unit length. These embeddings are particularly useful in the context of discrete 
holomorphy 3 ' 4 ' 5 , where the isoradiality condition allows the definition of discrete analytic 
functions. 

The intersection properties that we need to impose on the zig-zag paths of a planar 
lattice, to be able to define a Z-invariant Ising model, are equivalent 11 to the existence of 
isoradial embeddings of that lattice. Therefore a Z-invariant Ising model can defined on a 
planar lattice G if and only if G has isoradial embeddings. 

This result can better understood by considering a geometric interpretation of the model. 
In Section V we will show that a Z-invariant Ising model, with non-vanishing ferromagnetic 
couplings, on a planar lattice defines an isoradial embedding of the lattice which is character- 
ized by a set of geometric angles determined by the model's coupling constant at criticality. 
This geometric interpretation is closely related with Baxter's geometric solution 2 differing 
only in the fact that the angles involved characterize the embedding of the lattice and not 
the rapidity lines themselves, which need not to be embedded in the plane as straight lines. 

The reverse statement is also true, given an isoradial embedding of a lattice one can always 
define a Z-invariant Ising model on that lattice using the geometric angles characterizing 
the embedding and an elliptic function parameterization 2 of the coupling constants. The 
Ising model on isoradial embeddings considered in references 3 ' 6 is the restriction to the 
critical point of such a Z-invariant Ising model. Its critical point coupling constants Jf, can 
be expressed in terms of the respective edge length Lij on the isoradial embedding of the 
lattice, 

C0sh2 jc= A. (1) 

This relation between coupling constants and geometry has been found previously in 
the study of geometric properties of regular lattices at criticality, as effective angles of the 
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corner magnetization 12 and shape dependent modular parameters 13 ' 14 ' 15 ' 16 , suggesting that 
the geometric interpretation reflects an important aspect of the Ising model critical behavior. 
The paper is organized as follows: 

In Section II we introduce the rapidity lattice Q and the diamond lattice Gq associated 
with an arbitrary planar lattice G. In Sections III and IV the basic results of the Z-invariant 
Ising model are revisited in an embedding independent way. In Section V we discuss the 
geometrical interpretation of the model. Finally in Section VI our conclusions are presented. 

II. THE RAPIDITY LATTICE 

To define the rapidity lattice Q associated with a given planar lattice G it is convenient 
to first introduce the diamond lattice of G, denoted by Gq. The diamond lattice 3 has as 
vertices the vertices of G and the vertices of the dual lattice G* and has edges connecting 
each vertex of G to the vertices of G* corresponding to its neighboring faces, an example 
is shown in Fig. 3. The faces of G<> are quadrilateral and surround an edge of G and the 
corresponding dual edge in G*. 

A zig-zag path on the edges of G, turning alternately maximally left and right, corre- 
sponds to a path of alternating diagonals on a sequence of quadrilateral faces of G$ in which 
each face connects with its two neighboring faces by two opposite edges, see Fig. 3. To 
each zig-zag path we associate a rapidity line, defined as the path in the dual of G<> which 
is contained by the corresponding sequence of faces. In the example of Fig. 3 three such 
rapidity lines are shown. Each edge of G has two associated rapidity lines in one to one 
correspondence with the two zig-zag paths which intersect at that edge. 

The rapidity lattice Q associated with the planar lattice G is formed by all the rapidity 
lines defined in this way. In general Q is not embedded in the plane as a set of straight 
lines. The calligraphic type is a reminder that we are not interested in Q as lattice or as an 
embedding of a lattice in the plane but as a set of intersecting curves each with an associated 
free parameter, the respective rapidity. In section V a geometrical interpretation of this set 
will be given, we will see that under mild restrictions the intersection pattern of Q together 
with its rapidities is a complete description of an isoradial embedding of G in the plane. An 
isoradial embedding of G, in which all its vertices are at unit length from the neighboring 
vertices of G*, corresponds to a rhombic embedding of G$ in which all the lattice faces are 
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rhombi of unit edge length. 



III. Z-INVARIANCE 

Given an arbitrary planar lattice G, the associated rapidity lattice Q can always be 
constructed using the procedure of the previous section. However a Z-invariant Ising model 
can be defined on G only if the rapidity lines satisfy certain additional topological properties. 
In this section we review the Z-invariant Ising model in an embedding independent way 
suitable for the discussion of these properties. 

The Ising model on an arbitrary planar lattice is defined in the usual way, by assigning a 
spin <7j = ±1 to each vertex i and a coupling constant Jij to the edge connecting the vertices 
i and j. The model partition function is defined as 



H \<ij> J 



(2) 



The model is said to be Z-invariant 2 if a star-triangle deformation of the lattice, as shown 
in Fig. 4, changes the partition function only by a multiplicative constant Z star = R Z triang i e 
with R being determined by the local coupling constants. 

Z-invariance places stringent conditions on the model coupling constants which can how- 
ever be solved 17 by introducing a parameterization in terms of elliptic functions of modulus 
k with < k < 1. Following reference 2 the coupling constants are parameterized in terms 
of an edge dependent argument = Jijij) accordingly with table I, by choosing a regime 
and a modulus k. The corresponding Ising model is Z-invariant if the edge arguments r^- 
satisfy the following conditions throughout the lattice 

the sum of the edge arguments around a vertex is 2A, (3) 
the sum of the edge arguments in the boundary of a m-sided face is {m — 2)A 

with the multiplicative constant being, for the vertex labeling of Fig. 4, 

/ 2^ y 

\ sinh 2 sinh 2 Jj^ sinh 2 J^j / 

For a straight line rapidity lattice it has been shown 2 that the conditions (3) are satisfied 
for edge arguments determined by the difference of the rapidities on the two straight lines 
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crossing each edge. This result is expected not to depend on the particular embedding of 
the rapidity lines in the plane but only on their topological properties. We will show that 
the straight line requirement can be relaxed and that the only condition we need to impose 
on the rapidity lines is that they do not intersect each other more than once or self-intersect. 

Using oriented rapidity lines 19 the argument r^- at the edge is chosen, in an embedding 
independent way, according to the orientation of the edge relatively to the two rapidity lines 
crossing it. For un-oriented edges there are two possible orientations, shown in Fig. 5, and 
the corresponding edge arguments are chosen to be 



where a and (3 are the rapidities of the two rapidity lines involved. 

This choice of edge arguments satisfies the conditions (3) everywhere on the lattice pro- 
vided that the rapidity lines do not intersect each other more than once or self-intersect. 
This result can be argued as follows. 

As in reference 1 we assume that there a rapidity line ao which is crossed by all other 
rapidity lines on for % > 0. The rapidity lines are oriented is such a way that if we look 
along the orientation of a the rapidity lines a^ cross a from right to left. We assume also 
that the lattice G is formed by the intersection of the rapidity lines after they cross the 
reference line ao- These assumptions do not place any restriction on the lattice G. 

Consider a face of Q surrounding a vertex of either G or G* and let its boundary edges 
be oriented as the respective rapidity lines. For our choice of rapidity lines orientations it 
follows that all the boundary edges are oriented away from one of the boundary vertices and 
towards a second one, represented in the examples of Fig. 6 by the downward and upward 
triangle respectively. We will call such an orientation of the boundary edges standard. 

Non-standard orientations are not allowed by the restriction that the rapidity lines do 
not intersect each other more than once or self-intersect. Consider for instance a face of Q 
formed by the intersection of five rapidity lines, with the boundary edges oriented as shown 
in Fig. 7. Each rapidity line crosses the reference line ao, intersects two or more other 
rapidity lines and moves towards infinity or to the boundary of a compact region enclosing 
the lattice G. Given the first rapidity line a±, subsequent rapidity lines will have to intersect 
a to the left of ai to avoid double-intersections and self-intersections. This can be done for 




A + a — (3 for orientation b) 



for orientation a) 



(5) 
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all lines except the last one, a^, which should cross ao at some point P' or P" and then reach 
some distant point Q'. Neither of these two requirements can be met without crossing 
more than once ct\ or a 4 or self-intersecting and this non-standard orientation is not allowed. 
A similar argument can be made against orientations with more than one boundary vertex 
for which the boundary edges converge or diverge. 

The reader can verify that for the standard orientations of Q the edge arguments defined 
by (5) satisfy the conditions (3) everywhere in the lattice and the resulting Ising model is 
Z-invariant. If the rapidity lines, or equivalently the lattice zig-zag paths, are allowed to 
intersect each other more than once or self- intersect, then the conditions (3) will not in 
general be satisfied and a Z-invariant Ising model cannot be defined. 



IV. THERMODYNAMIC PROPERTIES 

Using the standard arguments of Z-invariance 1 ' 2 the partition function, for real positive 
coupling constants, can be evaluated in the thermodynamic limit of the lattice G. For 
completeness we reproduce here these results. The partition function is given by 2 



z = (4n)i Jv ex P (x;^)) 

\<ij> J 



(6) 



where N is the number of vertices on the lattice and 



j., s 1 r / x ( x K'q'{x)\ J 

with the various quantities being defined in table I for the three distinct regimes. 

The local magnetization can also be evaluated 2 and is found to be, for any spin on the 
lattice, 

((1-Q- 2 ) 1/8 forfi 2 >l 
<*>={ ■ (8) 

10 forfi 2 <l 

The regime I corresponds therefore to a low temperature ordered phase, while regimes II and 
III correspond to high temperature disordered phases. A phase transition occurs between 
regime I and II for ft — 1. 
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V. GEOMETRIC INTERPRETATION AT CRITIC ALITY 



The intersection properties that the rapidity lines, or equivalently the lattice zig-zag 
paths, need to satisfy in order to define a Z-invariant Ising model on a given planar lattice 
are equivalent 11 to the existence of isoradial embeddings of that lattice. In this section 
we will establish this relation between the Z-invariant model and isoradial embeddings by 
considering a geometric realization of the model. 

At criticality the elliptic modulus k is zero, the complete elliptic integral A becomes 
7r/2 and the coupling constants, see table I, are parameterized by elementary trigonometric 
functions 

sinh 2 Jy = t an . ( 9 ) 

For non-vanishing ferromagnetic coupling constants with edge arguments in the interval 
0<Tjj <7r/2 the conditions (3) can be written in terms of geometric half-angles r^- = 9^/2 
as 

% = 2ix around each vertex of G (10) 
(it — 9ij) = 2ir around each face of G 

where the sums are respectively over all edges which meet at a given vertex of G and over 
all edges in the boundary of a given face of G. 

A set of angles satisfying these conditions characterizes a rhombic embedding of G$, in 
which its quadrilateral faces are drawn in the plane as rhombi with unit edge length. The 
angle 0^- is the inner angle of the rhombic face which surrounds the edge [ij] , see Fig. 8 for 
an example. The equations (10) guarantee that the rhombi can be joined together, at the 
vertices of G and G*, in a consistent way without overlapping. 

In such a rhombic embedding the edges of G$ crossed by a given rapidity line are parallel 
and we can associated to each rapidity line an angle characterizing this direction. Any 
rhombic embedding of Go is in fact completely characterized 9 by the intersection pattern 
of Q, which determines the relative positions of the rhombi, and the direction of the edges 
crossed by each rapidity line. The condition that the rapidity lines do not intersect each 
other more than once, or self-intersect, is a consistency condition needed to ensure that all 
interior rhombic angles are positive. 
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The embedding of G associated with a rhombic embedding of C<> is a strictly convex 
isoradial embedding. Therefore if a Z-invariant Ising model can be defined on a given lattice 
that lattice admits isoradial embeddings in the plane. A set of critical edge parameters in 
the interval < < 7r/2 can always be obtained by assigning angles to the rapidity lines 
in the range ]0, tt[ which increase relatively to the order of intersection of the rapidity lines 
a, with the reference line cto- 

The reverse statement is also true, if a planar lattice G has an isoradial embedding in 
the plane then we can define a Z-invariant Ising model in G. This follows from the fact that 
zig-zag paths on an isoradial embedding do not intersect each other more than once or self- 
intersect. The corresponding Z-invariant Ising model can be defined with edge parameters 
determined by the rhombic angles of the embedding = XOij/ir. From (9) it follows that 
the model coupling constants at criticality can be expressed in terms of the respective edges 
lengths according to equation (1). The Z-invariant Ising model with such coupling constants 
is similar to the geometrical model pointed out by Baxter 1 ' 2 , the only distinction being that 
the angles associated to the rapidity lines characterizes not the rapidity lines themselves 
but the parallel edges of G§ crossed by them. The restriction to the critical point of this 
Z-invariant Ising model is precisely the Ising model considered by the authors of references 3 ' 6 
in the context of discrete holomorphy. 

VI. CONCLUSIONS 

In this paper we considered the Z-invariant Ising model on arbitrary planar lattices intro- 
duced without reference to a straight line rapidity lattice. We have shown that a Z-invariant 
model can be defined on such lattices if and only if the lattice satisfies the topological prop- 
erty that zig-zag paths do not self-intersect or intersect each other more than once. The 
associated rapidity lattice is then defined with rapidity lines in one to one correspondence 
with zig-zag paths. 

This topological characterization has the geometric counterpart of the existence of iso- 
radial embeddings of the lattice, establishing a close relation between the Z-invariant Ising 
model and the geometrical setting used in the mathematical literature to define discrete 
versions of analytic functions 3 ' 4 ' 5 . Our results can then be stated in a purely geometric way: 
a Z-invariant Ising model can be defined on an arbitrary planar lattice if and only if the 
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lattice admits isoradial embeddings in the plane. 

The Ising model on isoradial embeddings studied in references 3 ' 6 was shown to correspond 
to the restriction to the critical point of a particular Z-invariant Ising model. In this model, 
similar to Baxter's geometric solution, each rapidity line has an associated geometric angle 
characterizing the direction of the edges of diamond lattice G<> which the rapidity line 
intersects. The coupling constants of model can be expressed at criticality in terms of 
the inverse length of the respective edge on the isoradial embedding (1), reproducing the 
relation between critical coupling constants and geometry found in studies of the corner 
magnetization 12 and shape dependent modular parameters 13 ' 14 ' 15 ' 16 of regular lattices in 
non-trivial topologies. 
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FIG. 2: An arbitrary planar lattice G and the two zig-zag paths which intersect at the edge [ij]. 
Zig-zag paths on a planar lattice turn alternately maximally right and left and are said to intersect 
when share a common edge. 



n 

sinh2J(a) 
A 

q(a) 
r(a) 



regime I regime II regime III 



l/k' 
sn(a) 
cn(a) 

K 

01 (a) 
dn(a) 



k' 

/; / sn(a) 
cn(a) 

K 

9(a) 
dn(a) 



ik'/k 

ijj sn ( a ) 
dn(a) 

K-iK' 

9(a) 
cn(a) 



sn(a)cn(a) sn(a)cn(a) sn(a)dn(a) 
k 2 K' k 2 K' K> 



TABLE I: Parameterization of the coupling constants J and related quantities in terms of elliptic 
functions of modulus k, as given in reference 2 . The conjugate modulus k' = (1 — fc 2 ) 2 , the complete 
elliptic integrals K, K' and the elliptic functions sn(it), cn(u), dn(u), H(u) ,Hi(u), &(u) and @i(u) 
are defined in the usual way 18 . 
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FIG. 3: The lattice of Fig. 2 and the corresponding diamond graph G<> as defined on section II. 
The vertices of G are shown as circles and the vertices of the dual lattice G* are shown as squares. 
The edges of G$ are shown in doted line. The rapidity lines of Q are paths on the dual of G<> which 
are contained on a sequence of faces without turns. Three such rapidity lines are shown including 
the two associated with the represented zig-zag paths. 
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FIG. 5: The two possible relative orientations of an edge of G and its two rapidity lines and the 
corresponding edge arguments. The dotted line is the reference rapidity line. 




FIG. 6: Standard orientations of the rapidity lines on the faces of Q surrounding vertices of G and 
G*, shown using the notation of Fig. 3. The boundary edges are oriented away from a vertex in the 
boundary and towards a second one, denoted by the downward and upward triangles respectively. 
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FIG. 7: A face of Q with a non-standard orientation of the boundary edges. The rapidity curve 
«5 cannot be completed by joining points P with either P' or P" and Q with Q' without multiple 
intersections or self-intersection, which by assumption are not allowed. 
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FIG. 8: An isoradial embedding of a lattice G represented together with its dual G*, the diamond 
lattice G<> and the rapidity lattice Q. Each edge of G has an associated rhombic angle. The angle 
associated to each rapidity line characterizes the common direction of all edges of intersected 
by it. 
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